Positivity violation for the lattice Landau gluon propagator 
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We present explicit numerical evidence of reflection-positivity violation for the lattice Landau gluon 
propagator in three-dimensional pure SU(2) gauge theory. We use data obtained at very large 
lattice volumes (V = 80 3 , 140 3 ) and for three different lattice couplings in the scaling region (/3 = 
4.2, 5.0, 6.0). In particular, we observe a clear oscillatory pattern in the real-space propagator C(t). 
We also verify that the (real-space) data show good scaling in the range t £ [0, 3] fm and can be 
fitted using a Gribov-like form. The violation of positivity is in contradiction with a stable-particle 
interpretation of the associated field theory and may be viewed as a manifestation of confinement. 
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I. INTRODUCTION 

In recent years, there has been considerable interest in 
the possible violation of spectral positivity for QCD and 
in its relation to color confinement (for a review see |l(). 
Let us recall Q that the reconstruction of a Garding- 
Wightman relativistic quantum field theory from the cor- 
responding Euclidean Green functions is possible only if 
they obey the Osterwalder-Schrader axioms |[|. In par- 
ticular, the requirement of positive definiteness of the 
norm in Hilbert space is expressed in Euclidean space by 
the axiom of reflection positivity. For a generic 2-point 
function D (x — y), this axiom reads 



d A xd A v r(- 



■xo,x) D(x-y)f (yo,y) > 0, (1) 



where /(xo,x) is an arbitrary complex test function. 
The above condition implies the existence of a Kallen- 
Lehmann representation for D(x — y), which is neces- 
sary for interpreting the fields in terms of stable parti- 
cles. Thus, a violation of QJ implies that the Euclidean 
2-point function cannot represent the correlations of a 
physical particle. This may be viewed as an indication 
of confinement |l|. 

The relation between reflection positivity and Eu- 
clidean correlation functions can be made explicit by con- 
sidering the spectral representation P, |4| 



D(p) 



dm 



2 P( m2 ) 

p 2 + TO 2 



(2) 



for the Euclidean propagator in momentum space. Then, 
the statement of reflection positivity is equivalent to a 
positive spectral density p(m 2 ). This implies that the 
temporal correlator at zero spatial momentum D{t, p = 
0) can be written as 



C(t) = D(t,0) 



dio p(uj ) e 



(3) 



We note that for general spatial momentum p one would 
have u> = \/p 2 + m 2 . [In the particular case p = con- 



sidered here, the decay behavior of D(t, p) provides di- 
rect insight on mass-like properties associated with the 
fields.] Clearly, a positive density p(uj 2 ) implies that 



C(t) > 0. 



(4) 



Notice that having C(i) > for all t does not ensure the 
positivity of p(w 2 )- On the other hand, finding C(t) < 
for some t implies that p{uj 2 ) cannot be positive, suggest- 
ing confinement for the corresponding particle. 

For the gluon, the Landau propagator is predicted to 
vanish at zero momentum ja, la, |3, lj- III llOJ . This im- 
plies that the real-space propagator D(x — y) is positive 
and negative in equal measure, i.e. reflection positivity 
is maximally violated |a, |6J, |8j] . An infrared suppressed 
Landau gluon propagator has been obtained in several 
studies in momentum space [III C2- LLj H3 • Numerical 
indications of a negative real-space lattice Landau gluon 
propagator have been presented in the 3d SU(2) case 
[L2\. in the magnetic sector of the id SU(2) case at fi- 
nite temperature [l]| and, recently, in the 4d SU(3) case 
for one "exceptional" configuration 15]. In this work 
we verify this feature in detail (see Section m, using 
data obtained at very large lattice sizes for the 517(2) 
case in three space-time dimensions 14]. At the same 
time, we try to fit the numerical data in real space (see 
Section ITTT1) by considering a sum of Gribov-like propa- 
gators jfj, [jj. Let us recall that an excellent fit of the 
(momentum-space) propagator by a Gribov-like formula 
has been obtained for the equal-time three-dimensional 
transverse gluon propagator in Ad SU(2) Coulomb gauge 
u& and for the 3d SU(2) Landau case |14|. while in Ref. 
|4| the (real-space) transverse propagator has been fit- 
ted using a a Stingl-like formula in the Ad SU(3) Landau 
case. Also, several fitting forms have been considered in 
Ref. |17J for the gluon propagator in momentum space. 



II. VIOLATION OF REFLECTION POSITIVITY 

An explicit (numerical) proof of violation of the con- 
dition J2J ma y be difficult if the correlation function 



C(t) only becomes negative at relatively large values of 
t. In particular, this is the case if C(t) is of the form 
C(t) = e~ A * f(t) and f(t) is only negative for t > 1/A. 
In this case it is helpful to consider alternative quantities 
G(t) that are positive if p{u) 2 ) is positive. [Consequently, 
finding G(t) < implies violation of positivity for p(u> 2 ).] 
For example, one can define [18| the function 



G{t) 



dt 2 



lnC(i) 



C(t)C"(t)-[^(t)] s 
[C(t)]2 

Using Eq. © we can write G(t) as 0, Q 

G(t) = (( W - H) 2 ) , 



(5) 



(6) 



(7) 



where the averages denoted by ( ) are evaluated in the 
measure duj p(uj 2 ) e~^ '* . Clearly, if the density p{tu 2 ) is 
positive, so is G(t). Let us note that for C{t) = e~ xt f(t) 
one gets G(t) = d 2 \nf(t)/dt 2 , namely we get rid of the 
exponential factor e~ A * and it should be easy to check 
numerically if G(t) — or equivalently f(t) — is negative. 
The quantity G(t) could be of particular interest in a 
Ad study, since in this case it is more difficult to obtain 
good data for large time separations. In the case of a 
Gribov-like momentum-space propagator [?J 



D{p) = p 2 / (p 4 + M 4 ) 



one obtains the real-space propagator 5] 
C{t) = — / dpD(p)e- ipt 



e- Mt /^ [Mt t: 

cos — = + — 

2M VV2 4 



Then, using Eq. ©, it is easy to check that 



G(t) 



-M 2 



2 cos 2 



Mt 
V2 



(9) 

(10) 

(11) 



which is negative for all values of t. 

Notice that if C(t) is negative for some value of t we 
cannot evaluate its logarithm in Eq. (JSJ), while the ex- 
pression in Eq. 10 is always well defined for C(t) ^ 0. 

On the lattice, the real-space propagator can be eval- 
uated using 



N-l 



c(t) = -^J2 e ~ 27Tlk0t/ND ( k ^ 



(12) 



fco=0 



where N is the number of points per lattice side and 
D(k) is the propagator in momentum space. If the lattice 
action satisfies reflection positivity [19j , then we can write 
the spectral representation 



C{t) 



E-E„t 



(13) 
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FIG. 1: Real-space propagator C(t) as a function of i for 
coupling (3 — 5.0 and lattice volumes V — 80 3 (above) and 
V — 140 3 (below). Errors have been evaluated using boot- 
strap with 1000 samples. All quantities are in lattice units. 



where r n are positive-definite constants. Clearly, this 
implies that C(t) is non-negative for all values of t. 

As in the continuum, we can consider G(t) using Eq. 
© or, equivalently, the function G(t) [C(t)] 2 . This quan- 
tity can be easily discretized on the lattice by 



1 



G(t,a) = — [C(t)C(t + 2a) - C{t + af 



(14) 



where a is the lattice spacing. Indeed, in the continuum 
limit a -> , one obtains G(t, a) = G(t) C(t) 2 + 0(a 3 ). 
Furthermore, defining 



-E n t/2 



-E n {t+2a)/2 



(15) 



we can use the Schwartz inequality 18] to show that 
G(t, a) > for all values of t and a if the r n 's are positive. 
Also, considering the effective gluon mass 



m(t) 



log[C(t + a)/C(t)} 



(16) 
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FIG. 2: Plot of |C(t)| as a function of t for lattice volume 
V — 80 3 and coupling (3 = 6.0. For clarity, errors are not 
shown. All quantities are in lattice units. The solid line is 
drawn to guide the eye. 




FIG. 4: Scaling for the real-space propagator C(t) as a func- 
tion of t (in fm) for lattice volume V = 140 3 and couplings 
13 - 4.2 (x), 5.0(D), 6.0(O). Errors have been evaluated 
using bootstrap with 1000 samples. 
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FIG. 3: Plot of —G(t, 1) as a function of t for lattice volume 
V — 140 3 and coupling /3 = 4.2. Errors have been evalu- 
ated using bootstrap with 1000 samples. All quantities are in 
lattice units. 



one gets e ™(*) _ e m (*+ a ) = G(t,a) a 2 e m( - t+a ^ / C(t + a) 2 . 
Thus, if G{t, a) > one obtains m(t) > m(t + a), i.e. the 
effective mass should decrease when considering a larger 
time separation t. As discussed in |18J. an increasing 
effective gluon mass has been obtained already in the first 
numerical studies of the gluon propagator |20j , suggesting 
a violation of reflection positivity. Note, however, that 
Eq. i|16|) is ill-defined if C(t) changes sign. 

Here we use the 3d SU(2) Landau-gauge data pre- 
sented in [l4j in order to check if the conditions C(t) > 
and G(t, a) > are violated for the gluon propaga- 
tor. (The data have been analyzed using the bootstrap 
method with 1000 samples; we checked that results do 
not change when using 500 samples.) As explained in 



Ref. [IJ, we set the physical scale by considering 3d 
SU(2) lattice results for the string tension and the in- 
put value y/a = 0.44 GeV, which is a typical value for 
this quantity in the Ad SU(3) case. Since we consider 
h = c = 1, this implies 1 jvnT 1 = 0.4485^. 

We find that the real-space propagator C{t) is negative 
for several values of t, showing a clear oscillatory behavior 
(see Fig. QJ. In analogy with Ref. [l(j we also plot, in 
Fig- El the function |C(t)|: the spikes reveal the change 
of sign in the propagator C(t). Finally, in Fig. |31 we 
plot the function — G(t, 1): one can see that, as in the 
Gribov-like propagator, G(t, 1) is negative for all values 
of t. Thus, we find an explicit violation of positivity for 
the lattice Landau gluon propagator. Let us stress that 
this violation is clearly observable for the three lattice 
couplings and for the two lattice volumes considered. 



III. SCALING AND FITS FOR C(t) 

It is important to check if the behavior obtained for 
C(t) satisfies scaling for the lattice parameters consid- 
ered here. To this end, we apply to the data the match- 
ing procedure described in |14l Section III] and consider 
t in physical units using [lj, Table 2] . We obtain that all 
propagators become negative at t rs 0.7 fm and that the 
minimum is reached at t m i n « 1 fm (see Fig. QJ . More- 
over, finite-size effects seem to become important only at 
t > 3 fm. This means that our data for t G [0,3] fm 
are essentially infinite-volume continuum results. Note 
that the Gribov-like propagator C(t) in Eq. I|10|l has its 
minimum at t min = n/(M y/2). Thus, the above result 
for t m in would imply M w ir/y/2 fm^ 1 sa 2.22 fm~ x — 
438 MeV = 0.995^/0-. Let us also observe that the 
momentum-space Gribov-like propagator D(p) [see Eq. 
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FIG. 5: Fit of C(t) as a function of t (in fm) using a sum of 
two functions of the type 117t for lattice volume V — 140 3 and 
coupling /3 = 4.2. We also display Ci(t) and C'2(t) separately. 



©] has its maximum at p m ax = M. In Ref. [lj] we 
obtained p max = OStoiV^ = 3501^° MeV. 

As said above, in Ref. [lj] we have fitted the gluon 
propagator in momentum space using a modified Gribov- 
like or Stingl-like formula with four or five parameters. 
We now try to fit the data in real space using 



C(t) = ce- xt/ ^cos(b + Xt/V2\ 



(17) 



which is a generalization of the Gribov-like propaga- 
tor in Eq. (|1€)|1 . Clearly, this function also corresponds 
to a G(t) that is always negative. Let us stress that 
only for b — 7r/4 does this propagator correspond to 
the Gribov-like propagator in Eq. JHJl. In particular, 
for b = one gets the momentum-space propagator 
D{p) = (p 2 + M 2 ) / (p 4 + M 4 ) , which is finite at zero 
momentum. As reported in |l4|. it is still not clear from 
our data (on a 140 3 lattice) if the zero-momentum gluon 
propagator vanishes in the infinite- volume limit, as pre- 
dicted in dHH. 

We fit the data obtained for the two largest physical 
volumes, i.e. V = 140 3 and [1 = 4.2, 5.0, with t in the 
range [0, 3] fm. As can be seen in Fig. \5\ the data are 
well fitted using a sum of two functions of the type l|17|l . 
The corresponding fitting parameters are reported in Ta- 
ble [I] The averaged mass scales are Ai = 1.69(l) v / a : = 
745(5) MeV and A 2 = 0.74(1) Jd = 325(6) MeV. One 
can also obtain good fits of the data in the whole t range 
by considering the Fourier transform of the sum of three 
Stingl-like propagators in momentum space. These fits 
(using 12 parameters) have been reported elsewhere |21| . 
It is evident that fits of the gluon propagator in real space 



(see also U) require more parameters than fits in momen- 
tum space. This is due to the fact that the infrared data, 
for which the modeling is still not well understood, are 
spread over the whole time interval by the Fourier trans- 
form done in the evaluation of the temporal correlator 
C(t). 

Recently, it has been suggested |22l l23j that the vio- 
lation of spectral positivity in lattice Landau gauge be 
related to the quenched auxiliary fields used for gauge 
fixing. We note that the fitting form proposed for C(t) 
in [23| (also considering 5 fitting parameters) describes 
reasonably well our data up to t — 3 fm — yielding a 
light-mass scale of about LlA^/a — 500 MeV — but can- 
not account for the oscillatory behavior observed at very 
large separations. 



P C\ Al C2 


62 A2 


4.2 0.368(6) 3.83(4) 0.70(3) 


0.099(6) 1.54(4) 


5.0 0.361(6) 3.72(3) 0.56(3) 


0.089(6) 1.75(5) 



TABLE I: Fit of the data using a sum of two functions of 
the type fT7ll . setting b\ — 0. We obtain x/d.o.f. « 0.24 
(respec. 0.19) for /3 = 4.2 (respec. 5.0). The number of d.o.f. 
is 13 (respec. 17). The values of Ai and A2 are in /m _1 . 
The relatively small x/d.o.f. is probably due to the use of the 
diagonal part of the covariance matrix only. 



IV. CONCLUSIONS 

Using data from the largest lattice sides to date, we 
verify explicitly (in the 3d case) the violation of reflec- 
tion positivity for the SU{2) lattice Landau gluon prop- 
agator. This is one of the manifestations of confinement 
discussed in [l| . For very large separations (t > 3 fm) 
the propagator shows a clear oscillatory behavior, but of 
course one needs a careful extrapolation to infinite vol- 
ume in order to verify if this behavior survives in that 
limit. In the scaling region, the data are well described 
by a sum of Gribov-like formulas, with a light-mass scale 
M « 0.74 v / ct = 325 MeV, where a is the string tension. 
As a final comment, one should always bear in mind that 
the Gribov-like propagator may not represent the true 
analytic structure of the gluon propagator, but it is il- 
lustrative of a possible mechanism of confinement for the 
gluons (see also the discussion after Eq. (18b) in 24]). 
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